Polarized Hyperons from pA Scattering in the Gluon Saturation Regime 
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We study the production of transversely polarized A hyperons in high-energy collisions of protons 
with large nuclei. The large gluon density of the target at saturation provides an intrinsic semi-hard 
scale which should naturally allow for a weak-coupling QCD description of the process in terms of 
a convolution of the quark distribution of the proton with the elementary quark-nucleus scattering 
cross section (resummed to all twists) and a fragmentation function. In this case of transversely 
polarized A production we employ a so-called polarizing fragmentation function, which is an odd 
function of the transverse momentum of the A relative to the fragmenting quark. Due to this fct-odd 

ff") . nature, the resulting A polarization is essentially proportional to the derivative of the quark-nucleus 

cross section with respect to transverse momentum, which peaks near the saturation momentum 
scale. Such processes might therefore provide generic signatures for high parton density effects and 

s ! ■ for the approach to the "black-body" (unitarity) limit of hadronic scattering. 

a 



It has been known for over 25 years that A's produced in collisions of unpolarized hadrons exhibit polarization 
perpendicular to the production plane. As of yet, such data are not available for very high energies where one expects 
that hadronic cross sections are close to their geometrical values (the "black body limit"). However, the BNL-RHIC 
collider will soon collide protons and deuterons on gold nuclei at energies of ~ 200 GeV in the nucleon-nucleon center 
of mass frame; later on, much higher energies will be accessible at the CERN-LHC. In this letter, we demonstrate 
that the polarization of A hyperons produced in the forward region in high-energy collisions of protons and heavy 
nuclei may generically be a sensitive probe of high-density effects and gluon saturation in the target. 
. The wave function of a hadron (or nucleus) boosted to large rapidity exhibits a large number of gluons at small x, 
which is the fraction of the light-cone momentum carried by the gluon. The density of gluons is expected to saturate 
, when it becomes, parametrically, of the order of the inverse QCD coupling constant a s [1]. The parton density at 
saturation is denoted by Q 2 S , the so-called saturation momentum. This provides an intrinsic momentum scale [2] 
which grows with atomic number and with rapidity because more gluons can be radiated in the initial state when 
phase space is big. For sufficiently high energies and/or large nuclei, the saturation momentum Q s can become much 
I | larger than Aqcd, such that weak coupling methods are applicable. 

Forward A production in pA collisions is dominated by high- a; quarks from the proton traversing the high gluon 
density region of the heavy nucleus. The quarks typically experience interactions with momentum transfers of the 
order of the saturation momentum. Thus, for large gluon densities in the target, such that the saturation momentum 
is in the perturbative regime, Q s > 1 GeV, the coherence of the projectile is lost, and the scattered quarks (having an 
average transverse momentum proportional to Q s ) fragment independently [3]. While nonperturbative constituent- 
s' ■ quark and diquark scattering and hadronization models [4] have been employed to understand hyperon polarization in 
collisions of protons with dilute targets, we expect that in the high-energy limit the presence of the intrinsic semi-hard 
scale Q s should naturally allow for a weak-coupling QCD description of the process. One can thus calculate the 
cross section for qA scattering in this kinematical domain within pQCD [5], and the deflected, outgoing quark will 
subsequently fragment into hadrons, which is described by a fragmentation function. 

In order to explain the transverse A polarization in unpolarized hadron collisions within such a factorized pQCD 
description, it has been suggested that unpolarized quarks can fragment into transversely polarized hadrons, for 
instance A hyperons. The associated probability [6,7] is described by a so-called polarizing fragmentation function, 
sometimes also called Sivers (effect) fragmentation function. Its main properties are that it is an odd function of the 
transverse momentum relative to the quark, fc t , and that the A polarization is orthogonal to kt, because of parity 
invariance. The polarizing fragmentation function is defined as [7] 1 : 

& N D h y q (z,kt) = D h y g (z,kt) - D h i /q (z,k t ) = D h y q (z,k t ) - D h t /q (z,-k t ), (1) 

and denotes the difference between the densities D h r/ q (z, kt) and D h i/ q (z, kt) of spin- 1/2 hadrons h, with longitudinal 
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1 Another commonly used notation for the polarizing fragmentation function is D± T , but with a slightly different definition [6]. 
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momentum fraction z, transverse momentum kt and transverse polarization f or |, in a jet originating from the 
fragmentation of an unpolarized parton q. Clearly, this fct-odd function vanishes when integrated over transverse 
momentum and also when the transverse momentum and the transverse spin are parallel. In order to set the sign 
convention for the A polarization we define 

A N D hyq (z, kt) e, A N D ht/q (z, \k t \) (2) 

\qxk t \ 

where q is the momentum of the unpolarized quark that fragments and Ph is the direction of the polarization vector 
of the hadron h (the f direction) . Fig. 1 shows the kinematics of the process under consideration and indicates the 
direction of positive A polarization for each quadrant in the A production plane. 




FIG. 1. Kinematics of the pA^AX process. The direction of positive A polarization is indicated for each quadrant in the 
A production plane. 

It should be emphasized that such a nonzero probability difference A N D h i / q (z, kt) is allowed by both parity and 
time reversal invariance. Generally it is expected to occur due to final state interactions in the fragmentation process, 
where the direction of the transverse momentum yields an oriented orbital angular momentum compensated by the 
transverse spin of the final observed hadron. This polarizing fragmentation function is the analogue of the so-called 
Sivers effect for parton distribution functions [8] , which yields different probabilities of finding an unpolarized quark 
in a transversely polarized hadron, depending on the directions of the transverse spin of the hadron and the transverse 
momentum of the quark. The Sivers effect can lead to single spin asymmetries, for instance in p^ p — > ir X , a process 
for which such (large) asymmetries have been observed in several experiments. 

Recently, such a single spin asymmetry in e p 1 ' — > e' it X has been calculated in a one-gluon exchange model [9] . 
Shortly afterwards it was understood [10] as providing a model for the Sivers effect distribution function. A similar 
calculation has recently been performed by Metz [11] for the production of polarized spin- 1/2 hadrons in unpolarized 
scattering, which can be viewed as providing a model for the polarizing fragmentation function. Here we will not 
employ such a model calculation, but rather use a parametrization for the polarizing fragmentation functions obtained 
from a fit to data [7]. However, these model calculations do demonstrate that nonzero Sivers effect functions can arise 
in principle. 

Due to the fct-odd nature of the polarizing fragmentation function it is accompanied by a different part of the partonic 
cross section (essentially the first derivative w.r.t. k t ) compared to the ordinary, unpolarized A fragmentation function, 
which is Ait-even. The characteristics of the resulting A polarization will turn out to be rather different from presently 
available data for hadronic collisions at moderately high energies and with "dilute" targets. These data show a A 
polarization that increases approximately linearly as a function of the transverse momentum l t of the A, up to It ~ 1 
GeV/c, after which it becomes flat, up to the highest measured l t values: /( ~ 4 GeV/c. No indication of a decrease 
at these high l t values has been observed. Furthermore, the polarization increases with the longitudinal momentum 
fraction £ and is to a large extent ^fs independent. These features do not change with increasing A [12-14]. The only 
A dependence observed is a slight overall suppression of the A polarization for large A and higher energies. For Cu 
and Pb fixed targets, probed with a 400 GeV/c proton beam [13,14], the magnitude of the polarization is about 30 % 
lower than for light nuclei. This effect is usually attributed to secondary A production through ir~N interactions [14]. 
The slight suppression shows no evidence for a dependence on l t in the investigated range 0.9 < It < 2.6 GeV/c, albeit 
with rather low statistical accuracy. It is clear that this data on heavy nuclei is not in the kinematic region where 
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saturation is expected to play a dominant role and the main differences to the results presented below are that in the 
saturation regime the transverse A polarization will depend on the collision energy and no plateau region is expected. 

We shall now present our calculation of A polarization in the gluon saturation regime, following ref. [7] regarding 
the treatment of the polarizing fragmentation functions. 

As mentioned above, in the calculation of the qA cross section one is dealing with small coupling if the target nucleus 
is very dense; however, the well known leading-twist pQCD can not be used when the density of gluons is large. Rather, 
scattering amplitudes have to be resummed to all orders in a 2 times the density. When the target is probed at a scale 
< Q s , scattering cross sections approach the geometrical "black body" limit, while for momentum transfer far above 
Q s the target appears dilute and cross sections are approximately determined by the known leading-twist pQCD 
expressions. 

At high energies, and in the eikonal approximation, the transverse momentum distribution of quarks is essentially 
given by the correlation function of two Wilson lines V running along the light-cone at transverse separation r t (in 
the amplitude and its complex conjugate), 

) • (3) 



° qA = J ^f S{q+ - p+) \W c tT |/ d 2 z t e^[V(zt)-l] 



Here, P + is the large light-cone component of the momentum of the incident proton, and that of the incoming quark 
is p + = xP + (q + for the outgoing quark). The correlator of Wilson lines has to be evaluated in the background field of 
the target nucleus. A relatively simple closed expression can be obtained [5] in the "Color Glass Condensate" model 
of the small-x gluon distribution of the dense target [2] . In that model, the small- a; gluons are described as a classical 
non-abelian Yang-Mills field arising from a stochastic source of color charge on the light-cone which is averaged over 
with a Gaussian distribution. The quark q t distribution is then given by [5] 



1 = "5+ M o + 7^2 C(qt) , 



dq+d 2 q t d 2 b P+ \ P+ J (2tt) 2 

C(q t )= J d 2 r t e^ jexp -2Q 2 S J (1 - exp(ip t • r t )) 



2 cxp 



d 2 Pt 1 



1 ■ (4) 



(2n) 2 pf. 

This expression is valid to leading order in a s (tree level), but to all orders in Q s since it resums any number of 
scatterings of the impinging quark in the strong field of the nucleus. The saturation momentum Q s , as introduced in 
cq. (4), is related to \-> tne total color charge density squared (per unit area) from the nucleus integrated up to the 
rapidity y of the probe (i.e. the projectile quark), by 

Q 2 s =4n 2 a 2 s ^t-l X . (5) 

In the low-density limit, \ is related to the ordinary leading- twist gluon distribution function of the nucleus, see 
for example [15]. From BFKL evolution, Q 2 evolves as ~ exp(Ay) <~ x A , with the intercept A ~ 0.3 [16]. Thus, if 
Q 2 ~ 10 GeV 2 at the proton beam rapidity (i.e. x = 1) and for A ~ 200 targets [5], then Q s ~ 3 GeV at x = 0.6, 
decreasing to Q s ~ 2 GeV at x = 0.05; furthermore, assuming Q 2 <~ A 1 / 3, scaling, then at x = 0.6, Q s drops from 
3 GeV to 2 GeV when the atomic number A of the target decreases from 200 — > 20. It is clear therefore that in order 
to be sensitive to high-density effects, experimentally one should study high-energy p A collisions in the forward region 
(where the polarization is largest anyway, see below) and with large target nuclei, and then compare to pp collisions. 
Below, we shall focus on polarized A production in a relatively small rapidity interval in the forward region, and so 
take Q s as a constant of order 2 — 3 GeV. 

The integrals over p t in eq. (4) are cut off in the infrared by some cutoff A, which we assume is of order Aqcd- We 
denote the momentum of the produced A by I = zq + k, with k the transverse momentum relative to the fragmenting 
quark. Assuming parity conservation in the hadronization process, only the component of k in the production plane 
contributes to the polarization Va, therefore in order to simplify the kinematics we choose k y = as was done in Ref. 
[7]. For forward kinematics, q + 3> qt, one then finds zq t ~ It — k t . The polarized cross section is given by 

= / f f,/p(",Q 2 ) J ^A"D A , / ,( 2 ,Q 2 ,t,)rtC(») (7) 

1 

= / dxp q/p (x,Q 2 ) J ^LA N D A r /q (f,<? 2 ,^) C(q t ) , (8) 
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where £ = l z /P z ~ xz is the longitudinal momentum fraction carried by the A. We assume that A N D A i / q (z, k t ) is 
strongly peaked around an average lying in the production plane, such that [7] 



/ 



d%A N D A y q (z,Q^k t )F(k t )~A N D A r /q (z,Q 2 )[F(k°)-F(-k°)] . (9) 



Note that is a function of z, see below. Alternatively, one could consider Gaussian distributions over k t [17], though 
the above simplified treatment is sufficient for our purposes. 

Considering the unpolarized cross section, we can safely neglect fc t -smearing from hadronization, which is of order 
Aqcd, while the quarks are typically scattered to much larger transverse momenta, namely of order Q s . In our 
numerical results shown below we also include the contribution from anti-quarks and gluons to the unpolarized cross 
section, although this is a small correction for £ > 0.1. For the polarizing fragmentation functions, only contributions 
from u, d, and s quarks (the valence quarks of the A) are considered [7]. Thus, we obtain 

$1 dxxf g/p (x, Q 2 )A N D A y q U,QA \C (f (l t - - C (f (Z t + fc°)Y 
Vk{IuO = - ; J ~^7 \ x V — • (10) 



$\ dxxf q/p (x, Q 2 )D A / q (|,Q2) C (f/ t ) 



The factorization scale is chosen to be the saturation momentum of the dense nucleus, Q 2 — Q 2 . A parametrization 
for A N D A -i i q {z) in terms of the unpolarized fragmentation function D A / q (z) was given in ref. [7]. It was obtained by 
performing a fit to available pA — > A^A" data (for light nuclei only), where the transverse momentum l t was required 
to be larger than 1 GeV/c, in order to justify the application of a factorized expression and of pQCD for the partonic 
cross section. Although doubts have arisen about the applicability of pQCD in the kinematic region covered by the 
available data [17], the resulting functions do exhibit reasonable features. Here, we shall employ those functions as 
an Ansatz to investigate the dependence of the A polarization on the saturation momentum Q s , which turns out 
not to depend on the detailed parameterization of the polarizing fragmentation functions. Rather, it is the fc t -odd 
structure (and the fact that it is peaked around an average nonzero transverse momentum) that is responsible for the 
dependence. Of course, future parameterizations can be easily implemented. 
To be explicit, we use 

A N D Ayq (z, Q 2 ) EE N q z^ (1 - *)* , (11) 

where 

N u = N d = -28.13, N s = 57.53, c q = 11.64, d q = 1.23. (12) 

The average transverse momentum acquired in the fragmentation is parameterized as 

k° t = 0.66 z - 37 (1 - z) 50 GeV/c. (13) 

For the unpolarized fragmentation function D A / q in eq. (11) the parametrization of ref. [18] is to be used; strictly 
speaking, that parameterization holds for the fragmentation into A + A. However, in the forward region (£ > 0.1), 
one expects V A+A « V\- Furthermore, the parameterization of [18] assumes SU(3) symmetry: D A i u = D A / d — D A / S . 
However, the polarizing fragmentation functions A N D A -\/ q reduce the flavor symmetry to SU(2), since N u ^ ^ N s 
(N u = Nd was imposed in [7] to reduce the number of fit parameters). But even though A N D A i j s > \A N D A i / u ^\, 
the overall A polarization in the process under consideration is in fact dominated by the valence-like quarks of the 
proton, not by the strange quark. 

The polarizing fragmentation function describes the probability of an unpolarized quark to fragment into a trans- 
versely polarized A. Here, no difference is made as to whether the A is produced directly or as a secondary particle, 
for instance as a decay product of heavier hyperons like the S° or This second type is usually expected to have 

a depolarizing effect, which means that the degree of polarization is higher for the directly produced A's. The number 
of directly produced A's is estimated to be roughly 75 % of the total, such that the depolarizing effect could be on 
the order of 30 %. The polarizing fragmentation functions of ref. [7] thus effectively account for the depolarizing 
effect from decays, since they were obtained by a fit to data that does not discriminate between direct and decay 
contributions cither. 
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FIG. 2. Transverse momentum distribution of the transverse A polarization. Left: at fixed longitudinal momentum fraction 
£ = 0.5 and varying target saturation scale, Q s = 2, 3 GeV, respectively. Right: For Q s — 2 GeV and various £. 



A numerical evaluation of eq. (10) is shown in Fig. 2, using the CTEQ5L LO parton distribution functions for the 
proton [19]. Generically, one observes that Va is negative (due to the fact that u and d quarks dominate); it first 
increases with transverse momentum, then peaks at l t ~ Q s , and asymptotically approaches zero again. The fact 
that V\ peaks at l t ~ Q s has its origin in the fct-odd nature of the polarizing fragmention function: from eq. (10), 
V A corresponds to the difference of the cross sections taken with "intrinsic" transverse momentum k® parallel and 
antiparallel to the quark transverse momentum q t . Since is small, V a is essentially proportional to the derivative of 
da qA /d 2 q t , the differential quark-nucleus cross section, which varies most rapidly at q t ~ Q s (see also eq. (17) below). 
Consequently, \Va\ exhibits a maximum at such transverse momentum. This conclusion is independent of the details 
of the polarizing fragmentation functions; only the fct-odd nature and the fact that they are strongly peaked about 
an average matters. In contrast, A^-even distribution and fragmentation functions only probe the qA cross section 
itself but not its derivative with respect to qt ■ 

The behavior of "Pa^*) is qualitatively rather different when the quark cross section is taken at leading twist. In 
that case, not only is the magnitude of the polarization larger, but moreover V a in the forward region peaks about 
small transverse momentum < 1 GeV. This can be understood by noting that the derivative of the q A cross section 
at leading twist peaks in the infrared, contrary to eqs. (4,17). For a more quantitative evaluation of polarized A 
production in the "dilute regime" (hadronic collisions far below the unitarity limit, e.g. pp collisions at RHIC) we 
refer to Ref. [7]. 

To understand the behavior of eq. (10) in more detail, consider first large transverse momentum, q t 3> Q s - Here, the 
last two terms of eq. (4) can be dropped, since they contribute only via a 5{qt) term. At large transverse momentum 
the phase factor exp(iq t ■ r t ) in eq. (4) effectively restricts the integral over d 2 rt to the region r t < l/q t <C 1/Q S ; the 
first exponential can then be expanded order by order to generate the usual power series in 1/q 2 . The leading and 
subleading twists are (see also [20]) 



C(q t 
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'4 



4 Ql . q t 

--j-log- 



A 



O 



(14) 



This expression is valid to leading logarithmic accuracy. The first term corresponds to the perturbative one-gluon 
t-channel exchange contribution to qg — > qg scattering [20]. To leading order in /l t , the polarization given in eq. (10) 
thus becomes 



f*dxxf g/p (x,Q 2 )A N D A t /q (|,Q 2 ) 


x- 4 


1 + W k 




dxxf q / p (x,Q 2 )D A /q (liQ 2 , 


x- 4 


i + i^f iog^ 

7TX Z~ & A 





(15) 



It is known [7] that the polarization (for large It) is a higher-twist effect, i.e. it is suppressed by powers of the 
"intrinsic" transverse momentum at hadronization, k®, over the external momentum scale It- Eq. (15) shows that 
despite a partial cancellation the first power-suppressed correction to the quark-nucleus cross section (the subleading 
terms in the square brackets) enhance Va at large It, in agreement with the behavior at l t > 5 GeV in Fig. 2. 
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Regarding the scaling of Vk at the peak, consider the quark-nucleus cross section for q t <~ Q s 3> A. Again, the 
last two terms of eq. (4) can be dropped, while in the leading logarithmic approximation the argument of the first 
exponential reads 

-^ l0S ^A +0( ^ 2) - (16) 

The phase factor effectively cuts off the integral at r t ~ l/q t <~ l/Q s , and so l/r t A is large. We therefore replace 
l/r t — ► Q s in the argument of the above logarithm, since it is slowly varying and formally makes the expression 
well-behaved at large rt- The remaining integral leads to 

^/ \ 47T 2 / nq? \ 

c{qt) * q^q./a exp l-QfioioTA J • ^ 

This approximation reproduces the behavior of the full expression (4) about q t <~ Q s reasonably well. Expres- 
sions (16,17) are useful only when the cutoff A <C Q s , that is, when color neutrality is enforced on distance scales of 
order 1/A > 1/Q S . If, however, color neutrality in the target nucleus were to occur on distances of order 1/Q S [21] 
then A ~ Q s and one would have to go beyond the leading-logarithmic approximation. 
From eq. (17), Va is given by (to leading order in fcj/Zt) 

g I?°p(- a- Jg./A 

Thus, at the peak Pa scales approximately with l/{Q sS /\ogQ s / 'A), as indeed seen in Fig. 2. The strong dependence 
on the target gluon density, as parametrized by Q s , is rather different from leading- twist perturbation theory. 

As mentioned above, there is a related fc t -odd effect in processes with one transversely polarized hadron in the 
initial state (the Sivers effect), which can lead to asymmetries in p^ p — > it X [22], for example. Since at RHIC 
polarized proton beams are also available (for recent preliminary p^ p — > it X data from STAR, see Ref. [23] ) , one 
could investigate the process p^ A — > it X in the saturation regime. Similar signatures should arise in that process as 
for p A — > A^ X pointed out here. 

In summary, we have studied transverse A polarization in p A collisions at high energies and with dense targets. The 
resulting A polarization is quite different from that observed in p A and pp collisions to date, which presumably did not 
probe the saturation regime yet. To study the high-density limit, we have performed a weak coupling analysis of the 
hard q A scattering, determined by the saturation momentum Q s , and described the unpolarized quark fragmentation 
into a transversely polarized A hyperon by the so-called polarizing fragmentation functions. We observe that the A 
polarization peaks at transverse momentum <~ Q 8 , where it also scales approximately as l/(Q sA /logQ s /A) and hence 
is collision energy dependent. Moreover, no plateau region for larger transverse momenta is present. These features 
are independent of the details of the polarizing fragmentation functions, but rather occur due to their fct-odd nature. 
Similar effects are expected in the process of p^ A — > n X in the saturation region. Both processes can be studied, in 
principle, at the BNL-RHIC collider, and perhaps in the future at the CERN-LHC. 
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